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Abstract 

Lattice QCD simulations provide a promising way to disentangle different interpretations of hadronic 
resonances, which might be of particular relevance to understand the nature of the so-called XYZ par¬ 
ticles. Recent studies have shown that in addition to the well-established naive quark model picture, the 
axial-vector meson /i(1285) can also be understood as a dynamically generated state built upon the KK* 
interaction. In this work, we calculate the energy levels of the KK* system in the /i(1285) channel in finite 
volume using the chiral unitary approach. We propose to calculate the loop function in the dimensional reg¬ 
ularization scheme, which is equivalent to the hybrid approach adopted in previous studies. We also study 
the inverse problem of extracting the bound state information from synthetic lattice QCD data and comment 
on the difference between our approach and the Liischer method. 
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I. INTRODUCTION 


The/i (1285) is aaxial-vector state with mass m = 1281.9±0.5 MeV and 
width r = 24.2 ±1.1 MeV yj]. In the naive quark model, this state is assigned as a = ^Pi 

state. In recent years, however, it has been suggested to be a dynamically generated state made 
from the KK* interaction, together with its axial-vector counterparts yilSD. Such a picture has 
been extensively tested in the past decade ll4l-ll4ll. All these studies yield consistent results that 
the ground-state axial-vector mesons can be understood as dynamically generated states or at least 
contain large pseudoscalar meson-vector meson components. 

Lattice QCD (LQCD) simulations can be applied to study the properties of hadrons from first 
principles using quark and gluon degrees of freedom. Although studies of ground-state hadrons, 
which do not decay via strong interactions, have been well established and have turned out to be 
quite successful (see, e.g.. Refs. (IlsL Ilblll. studies of resonances are more challenging, since they 
do not correspond to discrete energy levels on the lattice, and considerable additional efforts are 
needed to extract physical information from LQCD simulations. The Liischer method is the de 


facto standard one in the case of single channel two-body elastic scattering 11171 


ISh . In this frame¬ 


work, the discrete energy levels obtained in LQCD simulations are related to the scattering phase 
shifts in infinite spacetirneQ In Ref. [26], the authors have developed a new effective approach to 
connect the LQCD discrete energy levels with the physical phase shifts (energies) by keeping the 
full relativistic two-body propagator, from which the Liischer formulation can be derived. This 
new approach has been applied to study finite volume effects in the meson-baryon interaction in 


the liilich model [|27h : the KD, rjDs interaction 


33 ; 


the pion-kaon 


scattering OOl. 


3111 : the 


DN, vrSc interaction ll32ll : the vrp interaction 111 ih : the tttt interaction yjj; and the KN interac¬ 
tion |34]. 

In the present work, we apply this approach to study the KK* interaction in the /i(1285) chan¬ 
nel. The /i(1285) is peculiar in the chiral unitary approach since it is made from the single channel 
KK* interaction and is located below the KK* threshold. As a result, it appears as a bound state 
in the dynamical picture. Its experimental width can be obtained from considering other coupled 
channels (see, e.g.. Ref. [Il3n i without affecting its nature being dominantly a KK* bound state. 
Inclusion of high-order kernels in the chiral unitary approach is found to have negligible effects 
^ Although in the present work we only need to tackle a single channel problem, it should be noted that the Liischer 


method has been generalized to the case of multichannel 


and rjK channels has recently been done in Refs. ||23 


scattering ll9 


2211 . A thorough study of the coupled ttK 


2411 and for the coupled tttt and KK channels in Ref. 02-511 . 


















































on this picture [Il2n . 


II. THEORETICAL ERAMEWORK 


A. Continuum 


In the ehiral unitary approaeh of Ref. lU], the /i(1285) is made of a single ehannel -^{\K*K) + 


\K*K)). The relevant l/-matrix is 


y(s) = 


e • e 


8/2 


-3) 


3s - (M" + + M'2 + m'2) - -(M" - 


f/2 /2'i 

m 


( 1 ) 


where / is the pseudosealar decay constant, s the invariant mass squared, e (e') stands for the 
polarization four-veetor of the ineoming (outgoing) K*. The masses M (M'), m {m') eorrespond 
to the initial (final) K* and K, respeetively. The potential V is unitarized via the following Bethe- 
Salpeter equation BO : 

T=ll + VG]-\-V)e-e', 


( 2 ) 


where G = G(1 + and q is given by 


q = —p\/[s — (M + m)2][s — (M — m)^]. 
2y s 


(3) 


The sealar loop funetion G has the following form: 

d^q 1 


G(\/i) = i 


{P — qY — M2 + ie — rn? + ie ’ 
with P the total incident momentum, whieh in the eenter-of-mass frame is (-y/s, 0, 0, 0). 


(4) 


The loop funetion G is divergent and needs to be regularized. This can be done either in the 
dimensional regularization seheme or in the eutoff seheme. In the former, the loop funetion reads 


Df 1 r / N 1 3^2 - m 2 + S m2 

+ -^[ln(s — (M^ — m?) + 2g-v/s) + ln(s + (M^ — m?) + 2q^/s) 
y/s 

— ln(—s + (M^ — rrY) + 2qy/s) — ln(—s — (M^ — wY) + 2qyfs)]}. 


(5) 


In our work, the regularization parameters are ehosen to be a(/i) = —1.85 and jj, = 900 MeV y]]. 
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B. Finite volume 


To study the /i(1285) meson in finite volume, one replaces T of Eq. Q by T, obtained using 
the same equation with the same potential and replacing the G-function in Eq.© by its counterpart 
defined in a finite box of size L. The function G in finite volurne, G, can be calculated again either 
in the dimensional regularization scheme, the cutoff scheme Obll . or a combination of both—the 


hybrid approach [12811 . To remove small unphysical discontinuities in the cutoff scheme, a smooth 
cutoff has been implemented in Ref. [26]. In the hybrid approach [^, an average of the results 


obtained with several sharp cutoffs is taken. This can save computational time when very large 
cutoff values are used. 

In principle in finite volume one mixes partial waves due to the cubic, rather than spherical, 
symmetry of the finite boxes chosen in the lattice simulations. The problem has been thoroughly 
studied in Ref. (IlSll and it is particularly relevant when one performs lattice simulations for parti¬ 


cles in a moving fra me 1122 , 
is also done in Ref. 


ne [22 

m\. I 


35 I-I 44 I 1 . The formulation for moving frames along the lines of Ref. [l26h 


In the present paper we only study systems with the two particles at rest 
interacting with S'-waves. We shall discuss the mixing in detail in Sec. IV, but we anticipate that 
for the levels that we consider in the inverse analysis, only the single channel with E=0 is relevant. 

In this work, we propose to calculate G in the dimensional regularization scheme. Introducing 
the so-called finite-volume correction, 6G, G can be written as: 


G = G^ + 6G, 


Eor the loop function of Eq. (|4|) , 5G has the following form [l45ll : 


( 6 ) 


5G = G{L) - G(oo) = 


dx5^l2{M‘^{s)), 


(7) 


where 


= {x^ — x)s -f + (1 — x)rn? — ie. 


( 8 ) 


Depending on the value of = \fP^, G needs to be treated differently. 
M + m, Sr(A4^(s)) can be written as a sum of the following three parts 1146[. 


n the case of y/s > 

d): 




( 9 ) 
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where the 2,3 defined as 


r 

9i 


r 

92 

r 

93 


E { |*ge ^ ' 

f*°° f 1 

Jo 27r2 \ [^ + A^2(-g)]r 
dr(M^(m^J) — r(x^ — x)(s 


1 r(x^ — x)(s — m^g) 

1 r(a;^ — a;)(s — 1 

[g2 + M^{m^gg)Y ^ [g^ + >12(^2j]r+i J ’ 

-ml)6r^YM\ml)). 


, ( 10 ) 

( 11 ) 

( 12 ) 


The separation seale needs to satisfy < M+m = M^* +mK- In the case of ^/s < M +m, 
5r{M.‘^) has a much simpler form 11451] : 

0 -l/ 2 -r/ / KA 2 '\i- 2 r 

^r{M\s)) = - ^ 3 ),^ ^ -5^(LV:^|n|)-3/2+^iT3/2-.(i^V:^|n|), (13) 

^ (^1 n^O 

where Kniz) is the modified Bessel function of the second kind, and 


00 00 00 

E E E (i-< 5 (N.o)). ( 14 ) 

n^O nx=—oony=—oonz=—oo 

with n = {tlx, ny,nY- It should be mentioned that the discrete summations in Eqs. (110113!) are 
only taken up to a certain number, |n|inax = where L and a are the lattice size and lattice 
spacing, respectively. Nowadays, most LQCD simulations adopt a L/a in the range of 16 ~ 32. 

In the hybrid approach, the finite volume effect is calculated in the following way: 


6G = lim 

Qmax ^CXD 


Qmax 


IS E 


|•q<qr, 


qi 


d?q 

(27r)= 


m 


where the function /(g) is 


m = 




2 a;(g)a;'(g) — (a;(^ + + ie 

with q = (n G Z^), + (p, u'{pj = pM‘^ + p, and E = y/s. 

In the Luscher method, the function I(q) of Eq. (fT^ is reduced to 11261] 

1 1 


/(g) = 


2E p'^ — q^ + ie' 


(15) 


(16) 


(17) 


where p = X^/'^{E'^,M'^,m?)/2E. 

In the present paper we are also treating the K* as a stable particle, while in fact it has a width 
of around 45 MeV. In an unquenched calculation if one uses interpolators of K*K one would 
reach the decay channels and one would have to deal with the three-body channels of The 
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formalism to deal with three body-systems in finite volume is also available in Refs. ll48l - l50ll . For 
two-body systems with one unstable partiele, one ean use a formalism in whieh the self-energy of 
the unstable partiele is diseretized in the moving frame 111 ill . We shall not do this here, although 
when more refined lattice calculations are available it would be interesting to tackle this problem. 
There are reasons not to do that at the present time. One of them is that many of the present 
lattice simulations use large pion masses where the decay channels would be blocked, but even 
there they can assess the existence of a bound state of KK* nature. The second reason is that in 
present lattice simulations, even using unquenched calculations, levels tied to channels that couple 
to certain quantum numbers do not show up unless explicit interpolators for this particular channel 
are explicitly used as interpolators. This was the case on the system looking for the ai(1260) 


resonance in Refs, oi 


520 and in the KD system in Ref. 11290 . where the levels associated to the 


coupled channel rjDs also did not show up in the simulation. The reason for this fact seems to 
be that the coupled channels not considered would show up in the time evolution at times where 
noise appears in the simulations, preventing any signal from being seen. The argument has stronger 
weight for the decay channels of resonances with a small width, like the present one with T = 24 
MeV. Obviously, there would be problems in the interpretation of the levels if these depend on the 
interpolators used, but the idea is to use interpolators with maximum overlap with the actual states, 
and there the effective field theories that we are using are of much help since they are telling the 
nature of the states under consideration. Then we suggest using interpolators that accommodate 
this structure, and in the present case these would be KK* interpolators. 


It is true that the consideration of t 


leads to changes in the spectrum yT 


le decay channels of the particles involved in a problem 


,y, 


530 and that to get the proper spectrum multihadron 


interpolators should be used ll54l] . but also, as mentioned in Ref. ll54l] , one can and must restrict 
oneself to lower energies if the interpolators accounting for the inelastic spectrum are not used. 
Concerning the present case we can use the analogy of this work, where we have KK* and the 


K* can decay to iCvr, and the case of Ref. ill ih . where one had pvr and the p could decay to vrvr. 


In spite of the large width of the p, the first level was very similar in the analysis with a stable p 
or a decaying p. The second level changed a bit more in both approaches, but it is reasonable to 
expect that with a smaller width of the K*, the differences would be much smaller. This, and other 
reasons that we will discus in Sec. IV concerning partial wave mixing, advise us to make use of 
only the first two levels that we shall discuss in the next section. 
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III. RESULTS AND DISCUSSIONS 


A. The energy levels 


The left panel of Fig. [T] shows the energy levels as functions of the cubic box size L obtained 
in the dimensional regularization scheme. For the sake of comparison, we show as well the en¬ 
ergy levels obtained in the hybrid method with q^ax = 4000 MeV. With the scale of Fig. [H the 
two curves are hardly distinguishable. However, as noticed in all previous works, there are some 
unphysical discontinuities in the hybrid approach, which disappear with an average of the results 
obtained with several sharp cutoffs or with a smooth cutoff [26]. This can be better appre¬ 
ciated from Fig. [2l which shows that the dimensional regularization scheme exhibits no sign of 
fluctuation, where small fluctuations can still be seen at a cutoff value of about 7000 MeV in the 


cutoff (hybrid) approacho In the following, unless otherwise noticed, we work with the dimen¬ 


sional regularization method. 




FIG. 1. Energy levels of the KK* system with G obtained from the dimensional regularization scheme 
in comparison with those obtained in the hybrid approach (left) and the Liischer approach (right) with 
Qmax = 4000 MeV. The lattice size L is given in units of where ttItt is the physical pion mass. 


The energy levels obtained in the Liischer approach are shown on the right panel of Fig. [T] as 
functions of the cubic box size L, in comparison with those obtained in the dimensional regular¬ 
ization scheme. It is clear that at least for the two lower energy levels, the Liischer results show 
stronger fluctuations than those of the dimensional regulmzation approach (also than those of the 
hybrid approach). Furthermore, it is shown in Ref. 113ill that the deduced phase shifts from the 

^ In the dimensional regularization scheme, for the sake of comparison, ^max has been related to |n|niax via ^max = 


\ fl\ 

L‘'|max‘ 
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FIG. 2. Finite-volume correetions, 6G, for the KK* system caleulated in the dimensional regularization 
scheme and the cutoff scheme with a sharp cutoff, (?max> at the energy E = 1250 MeV (left) and E = 1400 
MeV (right) and with L = 2.5/m7r. 


Liischer method can deviate by about 20 percent from the effective approach of Ref. 11281] at the 


energy region where the resonance dominates, at least for the ttK interaction in the K* channel 
(see Fig. 12 of Ref. [31]). 


As discussed in Ref. [|26h . the new terms incorporated in Ref. [|26l] with respect to the Liischer 


approach are exponentially suppressed and one would wonder whether other exponentially sup¬ 
pressed contributions from t and u channels, neglected in both approaches are not equally relevant. 
In this sense, explicit calculations of these effects done for mesons in the scalar sector 115511 . or the 
vector sector obil . show them to be negligible for lattice sizes bigger than L = 1.5 rn~^. 


The ai(1260) and 6i(1235) states have recently been studied in = 2 lattice QCD 11521] . 
where in addition to qq interpolators, meson-meson interpolators were also taken into account. 
Compared with the ai(1260) and 6i(1235), the /i(1285) is more suited to test the dynamical 
nature of the axial-vector mesons because of the following reasons. First, it is a single channel 
problem. Second, it is a bound state. Therefore it appears as a discrete energy level even in LQCD 
simulations. Third, it is built from the interaction of two strange mesons, which makes it less 
susceptible to chiral extrapolations. 

The ground-state pseudoscalar mesons and vector mesons have been studied in a number of 


rif = 2 + 1 LQCD simulations liSTUblh. Some of the gauge configurations are available on the 


International Lattice Data Grid, e.g., the PACS-CS configurations llSSh . which in principle makes 
a study of the /i(1285) straightforward. In TablejH we show the masses of the /i(1285) calculated 
in our framework, defined as the energies where T has a pole below the KK* threshold, with the 
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K, K* masses, and the lattice size L of the PACS-CS configuration 11581] (note, however, that these 


masses are calculated there with only qq interpolators). It is interesting to note that the /i(1285) 
remains as a bound state at these unphysical situations and the binding energy increases as the 
masses of its components increase. 


TABLE I. Masses, M, and binding energies, B, of the /i(1285) at unphysical quark masses and in finite 


volume. The K and K* masses are those obtained by the PACS-CS nj- = 2 + 1 simulations 1581] . In the 


last row, the numbers in the parentheses are the uncertainties coming from those of the K* and K added in 
quadrature. All the energies are in units of MeV while the lattice size L is in units of fm. 


Inputs 

Confl Conf2 Conf3 Conf4 Conf5 Conf6 

Physical 

nriK 

554(8) 594(9) 582(9) 635(9) 713(10) 789(11) 

495.0 

Mk* 

939(17) 984(16) 963(16) 1015(15) 1078(17) 1156(17) 

893.1 

L 

2.90(4) 2.90(4) 2.90(4) 2.90(4) 2.90(4) 2.90(4) 

oo 

M 

1367 1442 1412 1506 1635 1785 

1286 

B 

126(19) 136(18) 133(18) 144(17) 156(20) 160(20) 

102.1 


In Fig. [3l we show the mass of the /i(1285) as a function of the lattice size L at six different 
combinations of light and strange quark masses, corresponding to those of the PACS-CS configu¬ 
rations. It is clear that the results already approach their continuum limits at a lattice size of two 
to three times 1/m^. 

B. The inverse problem 

In this section we tackle the inverse problem of extracting an effective potential from discrete 
energy levels of LQCD. Close to the KK* threshold, one can assume a potential of the following 
form: 

V = a + b[s — {Mk* + rriKy]. (18) 

The two parameters a and b can be determined by fitting to the lattice energies. 

We assume that the first and second energy levels shown in Fig. [T] are ”LQCD” data. We 
take three energies from the first level and three more from the second one, and assign them an 
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FIG. 3. Mass of the /i(1285) as a function of the lattice size L at six combinations of light and strange 
quark masses, from lower to upper corresponding to Confl to Conf6 of TablelH respectively. 


error of 10 MeV. The eorresponding values of L are: 2.0 3.0 4.0 Performing a 

least-squares fit, we obtain a Xmm ~ 1 x 10“^ and the following two values: 

a = -157 ±29, 5= (-1.4±1.1) X 10"^MeV"^ (19) 


With the potential of Eq. (fTSl) . by solving the eorresponding Bethe-Salpeter equation, one finds 
a bound state at M = 1286 ± 37 MeV, whose eentral value eoineides with the original value 
we started with. It should be noted that although the bound state approaehes its eontinuum limit 
as L inereases, the potential approaeh has the advantage that it ean eonneet the LQCD energy 
levels at moderate L or small L with the binding energies in the eontinuum in a quite aeeurate and 
model-independent way (for a relevant and extensive diseussion, see, e.g., Ref. ll62ll f. 

Of eourse, for the ease at hand, one does not need to go through the inverse proeess to obtain 
the /i(1285) beeause it appears as a bound state. Nevertheless, this proeedure allows us to obtain 


an effeetive potential in a more or less mode 


-independent way. 


Following the approaeh of Refs. ll63l - l65ll . one ean quantify the relative eontributions of the 


meson-meson eomponent in the /i(1205) wave funetion. The eoupling eonstant of a resonanee to 
its eomponent ehannel ean be ealeulated as follows, 


= lim (s - so)T = lim 

S^SQ S^SQ V — (jr 


dv- 


dG 

ds 


( 20 ) 


s=so 


10 





















where sq is the pole position. From the above equation, one ean obtain the identity, 

,dG ,dV-^ , 

os OS 


( 21 ) 


The first term gives the eontribution of the eomposite component being dynamically generated, 
while the second term gives the rest (e.g., genuine qq or missing meson-meson channels). For the 
/i(1285), we find that —g'^^ = 0.50, which implies that the meson-meson component accounts 
for about half of the /i(1285) wave function. Given the fact that the /i(1285) is located about 100 
MeV below the KK* threshold, this value does not seem that small. 


IV. T»-WAVES FOR KK* WITH CHIRAL LAGRANGIANS 


In dealing with the KK* system with finite volumes one has to look at spin projections and 
partial wave mixing. There has been much work done along these lines recently. By using qq 
interpolators [l66l-l69n new methods and suitable interpolators have been developed to project on 
the desired spin states. More relevant to our problem, using the Liischer formalism for scattering 
of two particles, several papers have dealt with this problem. A detailed study for the case of 
0“, 1/2+ interacting particles is done in Ref.J^, which is generalized in Ref. [71] to moving 
frames. The case of 1/2+, 1/2+ interacting particles is studied in Ref.172] ^ and applied to the 


deuteron case in Ref. 117511 . A formal extension to the scattering of particles with arbitrary spin 
is done in Ref. [l76ll . For a first study of coupled-channel effects in LQCD simulations, see, e.g.. 


Refs. 1123, 


2M- 


In the present case, we are concerned about the scattering of 0“ and 1~ particles in the rest 
frame of the particles. If we only took into account S'-wave interaction between the 0“ and 1“ 
particles, it can be shown that on the cubic lattice and for total momentum P = 0 the S'-wave only 
mixes with the G-wave. However, in infinite volume, the /i(1285) can decay into a pair of 1“ 
and 0“ particles via the P-wave. In this case, mixing of L = 0 and L = 2 can occur. In order to 
assess the relevance of this component in the problem that we study we go back to the theory that 
generates the interaction of these particles using chiral dynamics. The chiral Lagrangian for this 
interaction is given in Ref. 117711 by 

1 


C 


WPP — 


4 /^ 


;TT{[V^,d,V,][P,d’^P]), 


( 22 ) 


^ See Refs. ll73Ll74ll for the generalized Liischer formula in multichannel meson(baryon)-baryon scattering formulated 
in both non-relativistic quantum mechanics and quantum field theory. 
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which leads to the potential 


(23) 


V ~ (pi+P3)(P2+P4)e%, 


where Pi,P 2 ,P 3 , Pa eorrespond to the two ineoming and two outgoing momenta in KK* —)■ KK*, 
and is the polarization of the veetor. It is elear that this potential has no D-waves. However, 


D-waves are automatiea^ 


ly gene rated where this Lagrangian is reinterpreted by means of the loeal 


hidden gauge approaeh [I78l - l80ll and is generated by the exehange of a light veetor meson (p meson 


for the KK* interaetion). In this ease one has the explieit p propagator and the potential beeomes 


^ ^ (Pi+P3)-(P2+P4)6% 

-{pi-nY-ml ’ ^ ^ 

whieh this time develops a D-wave. It is easy to see that the ratio of the -D-wave to S'-wave is 
(using I Pi — Pal < mp for the derivation) 


% _ 2p^ 1 

Vq ~ 3m2 ^1^2 +pV2 


(25) 


with p the CM momentum. 

Now we look at the energy levels of Fig. 1, that we have used for the simulation. We aetually 
took the first two levels for > 2 in the inverse analysis. Then we eonsider the levels 2, 3, 4 
that have energies in the eontinuum and we find the ratios for LmT^ = 2 


0.002 

for level 2, 

(26) 

0.079 

for level 3, 

(27) 

0.208 

for level 4. 

(28) 


The numbers would be further redueed by Clebsh-Gordan eoeffieients of L = 2 and S' = 1 to give 
J = 1, whieh are unity for L = 0. For bigger values of L, these energies are smaller and these 
ratios also deerease. For instanee for = 3 and level 4 we would find V 2 /V 0 ~ 0.059. 

The diseussions eondueted here are useful, beeause sinee we have only used the levels 1 and 2 
for > 2, then we always have a ratio of C 2 /H smaller than two per thousand, and we ean 
safely ignore the mixing. However, we also see that if we were to use the level 4 in our analysis 
and for values of < 2 we would have mixture of the order of 25% whieh would require us to 
explieitly eonsider the mixing for a proper interpretation of the results. 

The D-wave deeay of the /i (1285) eould in prineiple induee more eomplieated mixing patterns. 
For spin-0 and spin-1 seattering, whieh is the present ease, it ean be shown that the J = 1 H-wave 
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of the /i(1285) does not mix with any of the various P- and F-waves. However, it mixes with 
the J = 3 F-wave. Indeed, at energies mueh higher than eonsidered here, J = 3 resonanees have 
been found that deeay into KK*, sueh as the 03(1850) yj]. 

However, for the purpose of a rough error estimate, we ean assume that the J = 3 F-wave 
is of a similar size as the F-wave induced by Eq. (24), such that the uncertainties quoted in 
Eqs. (26)-(28) might be larger by a factor of 2. If lattice data become more precise, a coupled- 
channel calculation including the 5'-wave and the two F-waves will be necessary. At nonzero total 
momentum, which is not considered here, the mixing can become more complicated 11711] . 


V. SUMMARY 

We have studied the KK* interaction in the /i(1285) channel in finite volume with the chiral 
unitary approach. The relativistic loop function was calculated in the dimensional regularization 
scheme and compared with the hybrid approach developed previously. It was shown that although 
both approaches yield the same results if treated properly, the dimensional regularization scheme 
is numerically more stable. In addition, we found that the Eiischer method fluctuates more strongly 
with the variation of the cutoff, but agrees with the hybrid method qualitatively. 

In anticipation of future lattice QCD studies, we have calculated the position of the /i(1285) 
at six different combinations of light and strange quark masses as a function of the lattice size L. 
If confirmed, this could provide another test of the /i(1285) being a dynamically generated state. 
Indeed, the KK* meson-meson component is found to account for one half of its wave function. 
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